International  Journal  of  Electrical  and 
Electronics  Engineering  Research  (IJEEER) 

ISSN  (P):  2250-155X;  ISSN  (E):  2278-943X 
Vol.  8,  Issue  2,  Jun  2018, 13- 18 
©  TJPRC  Pvt.  Ltd. 

HIGH  DYNAMIC  PRECISION  CONTROL  SYSTEM 

Y.  A.  VERSHININ 

Coventry  University,  Coventry,  England 

ABSTRACT 

This  paper  is  dealing  with  the  design  of  a  high  precision  adaptive  control  system  for  multivariable  systems  or 
processes.  It  is  the  well  -known  fact  that  the  desired  behavior  of  local  subsystems  in  a  multivariable  system  is  difficult  to 
obtain  due  to  the  interconnections  between  channels.  The  method,  presented  in  this  paper,  allows  one  to  decouple 
channels  and  to  adjust  the  desirable  performance  of  the  overall  system  according  to  reference  models  of  local 
subsystems.  The  stability  of  the  overall  system  is  guaranteed  following  the  stability  criteria  based  on  the  Lyapunov 
function 
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INTRODUCTION 

Many  industrial  processes  and  systems  include  several  inputs  and  several  outputs.  A  mathematical  model 
of  such  a  processor  system  can  be  represented  in  the  form  of  a  MIMO  system  (Many-Inputs-Many-Outputs)  [1] 
and  [2].  For  example  a  chemical  batch  reactor  with  several  inputs  and  outputs.  Another  example  is  an  airplane, 
where  the  pitch,  roll  and  yaw  angles  can  be  controlled  by  a  pilot.  Yet  another  example  is  a  robot  manipulator 
where  the  motion  of  each  joint  can  be  controlled  independently.  Due  to  the  interconnection  between  channels  via 
the  system’s  dynamics,  the  response  of  the  system  may  not  follow  for  the  desiredtrajectory  of  the  applied  input. 
In  order  to  solve  this  problem,  decoupled  control  algorithms  are  suggested  to  use.  However,  parameters  of  real 
systems  can  change  due  to  the  change  of  environmental  conditions,  for  example,  due  to  the  change  of  the 
temperature  of  the  environment.  Therefore,  the  outputs  of  real  systems  will  also  change.  It  is  suggested  in  this 
paper  to  employ  adaptive  control  algorithms  in  order  to  provide  the  system’s  behavior,  according  to  the  reference 
model  for  each  channel.  It  is  admitted  that  the  dynamics  of  different  channels  can  also  be  different. 
Therefore,  the  relevant  reference  models  for  each  channel  can  be  employed  in  the  multivariable  system. 
The  stability  of  the  overall  multivariable  system  is  guaranteed  according  to  the  obtained  adaptive  control 
algorithms  based  on  the  Lyapunov  stability  criterion. 

THE  TASK  FORMULATION 

A  multivariable  plant  is  described  by  the  system  of  matrix  differential  equations  [1]: 

Ai  =  AnX1+A12X2  + ...  +  AlnX  n  +t/[, 

X2  =A2  lX1  +A22X2  +  ...+  A2nXn  +  U2, 

Xn  =  An\X  1  +  An2X2  +-+An  nX  n  +U„, 
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where 


X2  —  [.^21 ,  ^22’"'’ -^2v  ]  ’ 


Xn=\-Xnl,Xn  2,...,Xnk]  . 


Xi,X2,...,  X n  —  are  vectors  of  states  of  the  plant , 

Au  =  An{t),  Ap  =  Au{t),...,  Ann  =  Ann(t )  —  are  matrices  of  variable  paramerers  of  the  plant , 
Ui,Uiv..,Un  -are  control  vectors. 


It  is  required  to  find  such  feedback  laws  that  autonomous  motion  of  plant's  coordinates  at  each  channel  will  be 
provided  according  to  the  following  reference  model: 


v  m  —  A  m  v  m  i  n 
X\  -AuXl  +C/J, 

v  m  _  a  m  y  m  ±  r 
A  2  — /i22A2  _i"Lj2’ 


Y  rn  _  A  m  Y  m  n 
A  n  -  ^nn  A  n 


where 

ym  _\  m  m 

yVi  LA!pAi2v 

y  m  _  v  m  m 

yv2  —  La21’a22’* 

m 

**’  Al/77  J  ’ 

m  -\T 

. . . ,  x2v  ]  , 

ym  _  r  tn  m 

m  it 

...,xnki  . 

X"' ,  Xf ,...,  X"'  —are  vectors  of  states  of  the  reference  mode/, 

A";,A22,...,  A"'n  —  are  matrices  of  paramerers  of  the  reference  mod  el, 

Gl,G2,...,Gn  —are  inputs. 

Synthesis  of  the  Basic  Loop’s  Structure 

The  feedback  controls  can  be  chosen  as  follows: 

Ul=G1+KuXl+KuX2+...+  KuXn, 

U2  =G2+K21X,+K22X2  +...+K2nXn, 

un  =Gn+KnlXl+Kn2X2+...  +  KnnXn, 
where 

Kn  =  Kn(t),  Kl2  =  Kl2(t),...,  Knn  =Knn(t)-are  matrices  of  the  feedback  variable  gains. 


According  to  the  equations  (1)  and  (3),  the  closed-loop  system  can  be  represented  as  follows: 


Impact  Factor  (JCC):  7.9878 


NAAS  Rating:  3.19 
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(4) 


It  can  be  seen  from  (4)  that  the  system  "plant  +  controller"  will  be  dynamically  decoupled  on  each  channel  if  the 
following  conditions  are  achieved: 


Aij  +K,}  =  0,  (i,j  =  l,...n,  i  *  j) 


(5) 


In  order  to  guarantee  the  desired  motion  of  the  system  according  to  the  reference  trajectories  (2),  it  needs  to 
provide  the  following  relation: 


(6) 


Thus,  for  the  accomplishment  of  (5)  and  (6)  the  following  should  be  provided: 


a12,0  +  ^12,0  ~~  0, 
a12,l  +  ^12,1  =0’ 


(7) 


(f^i  o  +^ilo)  -  ^1,0’ 

<Au+W=cCi’ 


where 

«11,0  =fln,o(0,fliu  =am(0,-,a,„a.-i=flmi-i(0-a«  variable  scalar  coefficients  of  the  plant  (1), 
a'iiO’aui’—’annk-i~are  constant  scalar  coefficients  of  the  mod  el  (2), 

k  no  =  kt  ,0  (0,  kt  ,|  =  kUA(t),...,k,mjc_1  =knnk_l(f)-are  variable  scalar  gains  of  the  controller  (3). 

Synthesis  of  the  Adaptive  Control  Algorithms 

The  variable  parameters  of  the  plant  (1)  can  be  represented  as  follows: 


An  =  Afi  +  AA,  i , 
A12  =  A12  +A^12’ 


Where 
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A.j  =  const,  (i,  j  =  l,...n),  A°  —are  the  constant  matrices  whose  coefficients  correspond  to  the 
nominal  operation  mod  e  of  the  plant, 

A Ajj  =  A Ajj  (t),  A Aij  —  are  deviation  of  the  plant  parameters  from  their  no  min  al  operation  mod  e. 
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The  controller  parameters  can  be  represented  analogously: 

Kn=K°n+AKn, 

Kl2  =  K°2  +  AKl2, 


K .  ~  A'"  I  :\K. 


where 


K °  =  const,  (i,  j  =  l....n),  K°  -  are  matrices  of  cons  tan  7  coefficients  of  the  controller, 

AKjj  =  AKfft),  A Ktj  -  are  matrices  of  adjusting  coefficients  of  the  controller. 

Thus,  the  motion  of  plant  coordinates  on  a  non-nominal  operation  mode  can  be  represented  as  follows: 


=XK+AA!y(f))X,.+t/,  (1  = 

j= i 


The  feedback  control  is: 


(8) 


Ui=Gi+Y,(K(!j+^Kij^Xj,  (*=1, (9) 

i= i 

Substituting  the  equation  (9)  into  (8)  we  can  obtain  the  equation  of  motion  of  the  system  "the  plant  +  the 
controller"  on  each  channel: 


x  i  =  X 1 C <  +  A Ay  (t))X  J  +  £  (K  °  +  AKy  (t))X  J  +  G, . 
j= 1  7=1 


The  model  motion  equation  on  each  channel  can  be  chosen  as  follows: 


X,m  =A”'x;"+£(A°+K°)XJ+Gi. 

7=1 

j*> 


Comparing  (10)  and  (11)  and  taking  into  account  that  Ajj  +  K’u  =  A'f ,  we  can  obtain: 


(10) 


(11) 


E,  =  A™ Et  +  (A A,  (t)  +  A K„  (t))X ,  +  (£  A Ay  ( t )  +  £  AKy  (t))X  J 


7=1 

j*i 


7=1 

j*i 


where 


E,=Xi-X? 


(12) 


(13) 


Impact  Factor  (JCC):  7.9878 


NAAS  Rating:  3.19 


High  Dynamic  Precision  Control  System 


17 


Let  us  introduce  the  following  notation: 


Yll=Yll(t)  =  AAu(t)  +  AKll(t), 

Yy  =  Yij  (?)  =  A Ajj  (f)  +  A K,J  (0,  (i  *  ;). 


(14) 


Thus,  the  error  equation  (12)  can  be  represented  as: 


n 


(15) 


7=1 

j*i 


Algorithms  of  adaptation  can  be  specified  in  the  following  form: 

—  AK  =  'Y  (16) 

dt 

Admit,  that  parameters'  deviation  AA  =  A A(t)  is  differentiable  in  time: 

A A(t)=R(t)  (17) 

dt 

Therefore,  according  to  (13)  -  (17)  we  obtain: 

E  =  AmE  +  YX,  (lg) 

Y  =  '¥  +  R. 

where  R  =  R(t ) 

The  equations  (18)  with  eq.  (2)  Describe  the  dynamic  motion  of  the  adaptive  control  system  with  the  reference 


model. 


Admit,  that  the  matrix  T*  is  a  function  of  error  deviation  E  and  time  t,  and  yY(E,t)  =0  at  t  =  0. 


In  this  case  the  state  and  parametric  error  motion 


E  =  0,  Y  =  0 


(19) 


With  the  admission  that  R(t)  =  0  according  to  the  hypothesis  of  quasi-stability,  is  the  solution  of  the  system  (18). 


It  is  suggested  to  use  the  second  method  of  Lyapunov  in  order  to  obtain  the  adaptation  algorithms  VT  from  the 
condition  of  stability  of  the  zero  solution  (19)  of  the  system  (18))  [3],  [4],  [5]  and  [6]. 


The  quadratic  function  V  can  be  chosen  as  follows: 


V  =  }Et  PE  +  tr(YYr), 


(20) 


where 


y  =  const 

P  —  is  simmetric  matrix. 
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Q  —  is  negative  definite  matrix 

A  time  derivative  for  the  function  V  is  obtained  as  follows: 

V  =  yET  PE  +  yET  PE  +  tr(YYT  +YYT)=  y(A"'E  +  YX)T  PE  +  yETP(A'nE+  YX)  +  tr(YYT  +YYT )  = 
y[{A"'E)T  +  (YX)t]PE+  yETP{A"‘E  +  YX)  +  tr(YYT  +YYT)  = 

y[Er  ( Am  )T  +XtYt]PE+  yETPA’"E  +  yEr  PYX +tr{YYT  +YYT)  = 
yET  (Am)r  PE  +  yXTYTPE+  yETPA"'E+  yEr  PYX  +  tr(YYT  +  YYT )  = 
yET((Am)T  P+PA"')E+  yXTYT  PE  +  yEr  PYX  +  tr(YYT  +  YYT) 

After  some  intermediate  matrix  manipulations  we  obtain  that: 

¥  =  -yPEX  7  (22) 

The  derivative  of  the  Lyapunov  function  is  represented  as: 

V  =  ]EtQE  (23) 

Taking  into  account  that  A7  is  Hurvitch  matrix,  we  obtain: 

V  <  0  (24) 

Therefore,  according  to  (18),  (16)  and  using  conditions  of  quasi-stationarity  of  the  system,  it  is  straightforward  to 
obtain  the  adaptation  algorithms  as: 

AK=-]PEXT  (25) 

It  follows  from  the  above  that  the  motion  of  the  system  (16)  with  the  algorithms  (25)  is  stable  and  conditions  (5) 
and  (6)  are  satisfied. 

CONCLUSIONS 

The  developed  in  this  paper  control  algorithms  allow  one  to  improve  the  performance  of  a  multivariable  process 
or  system  due  to  the  decoupling  of  channels  of  the  process  or  system  from  the  interaction  via  the  process’s  or  system’s 
dynamics.  The  desirable  performance  of  the  process  or  system  can  be  obtained  using  the  relevant  reference  models  for 
each  channel.  Adaptive  control  algorithms  allow  one  to  adjust  the  overall  system’s  response  in  the  case  of  parameters 
change.  The  stability  of  the  overall  system  is  guaranteed  according  to  the  Lyapunov  stability  criterion. 
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